Surface Roughness and Hydrodynamic Boundary Conditions 
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We report results of investigations of a high-speed drainage of thin aqueous films squeezed between 
randomly nanorough surfaces. A significant decrease in hydrodynamic resistance force as compared 
with predicted by Taylor's equation is observed. However, this reduction in force does not represents 
the slippage. The measured force is exactly the same as that between equivalent smooth surfaces 
obeying no-slip boundary conditions, but located at the intermediate position between peaks and 
valleys of asperities. The shift in hydrodynamic thickness is shown to be independent on the 
separation and/or shear rate. Our results disagree with previous literature data reporting very large 
and shear-dependent boundary slip for similar systems. 
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It has been recently well-recognized that the classical 
no-slip boundary condition 1], which has been applied 
for more than a hundred years to model macroscopic ex- 
periments, is often not applicable at the submicro- and 
especially nanoscale. Although the no-slip assumption 
seems to be valid for molecularly smooth hydrophilic sur- 
faces down to a contact it is now clear that 
this is not so for the majority of other systems. The 
changes in hydrodynamic behavior are caused by an im- 
pact of interfacial phenomena, first of all hydrophobicity 
and roughness, on the flow. A corollary from this is that 
a theoretical description based on the no-slip condition 
has to be corrected even for simple liquids. What, how- 
ever, still remains a subject of hot debates is that how 
to correct the flow near the interface, and what would be 
the amplitude of these corrections. 

The simplest and the most popular way to model the 
flow is to use a slip-flow approximation 0, which as- 
sumes that the slip velocity at the solid wall is pro- 
portional to the shear stress, and the proportionality 
constant is the so-called slip length. Such an assump- 
tion was justified theoretically for smooth hydrophobic 
surfaces 0, 0, E| and was confirmed by recent surface 
force apparatus (SFA) |^ and atomic force micro- 
scope (AFM) (5i dynamic force experiments. Despite 
some remaining controversies in the data and amount 
of slip (cf. 01), a concept of a hydrophobic slippage is 
now widely accepted. 

For rough surfaces a situation is much less clear both 
on the theoretical and experimental sides. One point 
of view is that roughness decreases the degree of slip- 
page [3jEiH|i unless the surface is highly hydrophobic, 
so that trapped nanobubbles are formed to accelerate the 
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FIG. 1: Schematic of the dynamic AFM force experiment. 

flow 0, . An opposite conclusion is that roughness 
generates extremely large slip [l7| . 

We believe our letter entirely clarifies the situation 
with flow past rough surfaces, highlights reasons for ex- 
isting controversies, and resolves apparent paradoxes. 

We use a s peci ally designed home-made AFM-related 
setup 0, l2(j| to perform dynamic force experiments 
on a nanoscale (Fig.^. Glass spheres of radius R ~ 20±2 
/im (Duke Scientific, PaloAlto, CA) were attached with 
UV glue on the top of the rectangular tipless cantilever 
(length L — 450 /im, width w = 52 //m, spring con- 
stants k ^ 0.22 N/m). The spheres were then coated 
with a gold layer (50 nm) using a layer of chromiun (3 
nm) to promote adhesion. For the planar substrate we 
used a silicon wafer coated with a gold layer (50 nm). 
Both surfaces were then treated with a 1 mM solution of 
11-amino-undecane thiol (SII-(CIl2)ii-NH2) for 24 hrs 
to produce a chemically bound SAM. The advancing 
and receding water contact angle on the thiolated planar 
surfaces were measured with a commercial setup (Data 
Physics, Germany), and were found to be (69 ± 1)° and 
(63 ± 1)°, respectively. These values are close to those 
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for surfaces used in |2l|. Imaging of thiolated surfaces 
with a regular AFM tip revealed the root-mean-square 
roughness over a 1/im x 1/im is in the range 10-11 nm 
for a sphere and 0.5-0.8 nm for a plane. The maximum 
peak-to-valley difference is less than 45 nm for a sphere 
(Fig. |2Jl and less than 2.5 nm for a plane. This (smooth 
against rough) geometry of configuration allows one to 
avoid large scatter in data at separations of the order 
of or smaller than the roughness size. Such a scatter 
would be inavoidablc for two rough surfaces (depending 
whether the rough sphere is falling on a tip or in a valley 
of a rough plane) . Cantilevers were then fixed in a holder 
with the variable tilt angle, and the intermediate posi- 
tion with the angle 9 = (10 ± 2)° was chosen. A planar 
substrate was placed onto the bottom of a Teflon cuvette, 
which was filled with ImM NaCl (99.99%, Aldrich) aque- 
ous solutions. Water for solutions was prepared using a 
commercial Milli-Q system containing ion-exchange and 
charcoal stages. The deionized water had a conductivity 
less than 0.1 x 10"^ S/m and was filtered at 0.22 /im. 
To measure force-versus-position curves the cuvette was 
moved towards the particle with a 12 /im range piezo- 
electric translator (Physik Instrumente, Germany). This 
translator was equipped with integrated capacitance po- 
sition sensors, which allows to avoid any creep and hys- 
teresis and provided the position with an accuracy of 1 
nm in closed-loop operation. During the movement the 
deflection of the cantilever was measured with an opti- 
cal lever technique. Therefore the light of a laser diode 
(1.5 mW, 670 nm) was focused onto the back of the gold 
coated cantilever. After reflection by a mirror, the po- 
sition of the reflected laser spot was measured with a 
position sensitive device (United Detectors, UK, active 
area 30 x 30 mm^). The total force was calculated by 
multiplying instantaneous cantilever deflection with the 
spring constant. The distance h between surfaces was 
calculated by adding the piezo displacement to the de- 
flection of the cantilever, so that h — corresponds to 
the contact [of tips of sphere's asperities with a plane] 
(see Fig. nj. We stress, that since our plane is smooth, 
we have no ambiguity in determining this zero of separa- 
tion. 

The AFM force balance incorporates both (concen- 
trated) force on the sphere and the drag on the can- 
tilever (Fig. P), so that the total force measured in 
the AFM dynamic experiment is 
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where Fg and are the surface and the hydrodynamic 
forces, correspondingly, acting on a sphere, and Fc is 
the force due to distributed hydrodynamic drag on the 
cantilever. 

The total force Ft is proportional to the instantaneous 
deflection of the end of the spring from its equilibrium 
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FIG. 2: An AFM image of an apex of a gold-coated sphere 
(top) and a schematic representation of rough sphere-plane 
interaction near the point of a contact (bottom). 

position multiplied by the spring constant k 0] 

Ft^k{h-{ho + vt)), (2) 

where ho is the initial separation between surfaces, and v 
is the driving speed of piezo (negative speed corresponds 
to approach). 

The hydrodynamic force F^ between a sphere and a 
plane can be written as 22] 



67r/xi?2 dh 
^'^ 



(3) 



where fi is the dynamic viscosity, —dh/dt is the relative 
velocity of the surfaces. To finalize the description of Fh 
we should define the expressions for a correction function 
/*, which will be discussed later. Note that when /* = 1, 
Eq.|31transforms to famous Taylor's formula (which, how- 
ever, never appeared in any of G.I. Taylor's publications 
as discussed in H^). 

The drag force on a cantilever Fc is given by Q 
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where 7 — {2R+h)/{Lsm6). Here i? is a constant, which 
reflects the contribution from the Stokes flow to the can- 
tilever deflection, and represents an adjustable (dimen- 
sionless) parameter to the model. 
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FIG. 3: Hydrodynamic force acting on a sphere at high speed 
of approach. Only each 2nd point is shown. From bottom to 
top are data (symbols) obtained at the driving speed -6 fim/s, 
-10 nm/s, -20 /im/s, -30 /im/s, and -40 fim/s. Solid curves 
show the calculation results obtained for the same speed, by 
assuming f* — 1. 



The surface force Fg is assumed to be unaffected by the 
relative motion of surfaces, and was taken to be the equi- 
librium force being a function of only h. It was obtained 
from low speed (below 1 /im/s) force measurements. At 
distances larger than 20-25 nm no interaction was de- 
tected. In other words, no electrostatic contribution can 
be seen despite relatively large (9.6 nm) Debye length 
in 1 mM NaCl solution. This suggests that the surfaces 
are uncharged. Similar observations have been made be- 
fore for some other classes of thiols |24[. The range of 
the jump distance was always 15 ± 2 nm. The contri- 
bution of contact deformation to the jump distance was 
estimated using the experimental values for the pull-off 
force 0.70±0.05 mN/m (with the correction to the hydro- 
dynamic interaction) and the values of Young's modulus 
of the UV glue (3 GPa), as it was the softest material 
in our system. We also ignore the possible plastic flat- 
tening of the gold (50 GPa) asperities. For these pa- 
rameters the central displacement due to a contact de- 
formation is of the order of 0.1 nm, so that it can safely 
be ignored. We have fitted the experimental results by 
assuming Fg = —AR/Qh?, taking the Hamaker constant, 
A, as an adjusting parameter. The value A — 5.4 x lO"^*' 
J is obtained from fitting and was further used in all cal- 
culations. 

Fig. 13 shows the hydrodynamics resistance force cal- 
culated by subtracting Fs and Fc from the total force 
measured. Theoretical curves obtained by a numerical 
solution of differential Eq.(^ in the assumption /* = 1 
are also included. Note that cantilever contribution was 
found to be neither small or negligible. The adjusting 
parameter B reflecting Stokes-like flow on a cantilever 
was about 30-34 for our geometry of configuration and 
slightly varied from one experiment to another. Mea- 



FIG. 4: Hydrodynamic force versus inverse of separation ob- 
tained at driving speed -30 /im/s. Symbols represent exper- 
imental data. Solid curves show the theoretical results ob- 
tained within Taylor's formula. Upper line shows the cal- 
culations results in the assumption of a constant speed of 
approach. Lower curve corresponds to the real approach- 
ing speed. Dashed curve is the calculation results within the 
model described by Eq.|H]with bs = 48 nm, and dotted curve 
- by Eq. |7|with rs = bs = 48 nm. The approximation Eq. |7| 
with Ts = 45 nm provides the best fit to the data. 



sured Fh is much smaller than predicted by Taylor's the- 
ory. The deviations from theory are clearly seen at dis- 
tances 100-200 nm depending on the driving speed. One 
can conclude therefore that there is clear impact of rough- 
ness on the film drainage. 

In Fig. ^ the hydrodynamic force measured at driv- 
ing speed — 30/im/s is compared with Taylor's equation 
(solid curves). The force is plotted versus the inverse of 
the surface separation. Note that even Taylor's equation 
does not result in a linear plot. This would be expected 
only at a constant approach velocity, which is not the 
case due to forces acting on the sphere. However these, 
caused by decrease of dh/dt, deviations from linearity are 
much smaller than required to fit experimental data. The 
general analytical solution for rough interfaces does not 
exist in the literature probably since such a problem is 
outside the scope of a lubrication approximation. Below 
we analyze some approximate models. 

We have calculated the force expected between smooth 
slippery surfaces. The slip length is assumed to be 
roughly equal to zero for a sphere and to the height of 
asperities bg (as defined in Fig. [2Jl for a sphere, which 
is close to its definition. In this case the correction for 
slippage takes the form 22] 
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An improved fit to the data is obtained when slip is per- 
mitted (dashed curve in Fig. 21 ■ For this speed of ap- 
proach the best fit to Eq. was found with bg = 48 nm. 
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Eq. provides an excellent fit at distances larger than 
~ 100 nm. However, since data still deviate from theo- 
retical predictions at smaller separations, it is clear that 
slippage does not necessarily represent the roughness. 

Therefore, we further assumed that stick boundary 
conditions remain valid, but are applied to a surface de- 
fined at a distance towards a sphere: 



(7) 



By adjusting the value of we found that the hypoth- 
esis of a shift in separation gives a perfect coincidence 
between data and theoretical predictions (dash-dotted 
curve in Fig. 2) obtained with = 45 nm). We remark 
and stress that Eq. 13 provides an excellent description of 
the data even when h is much smaller than . [Another 
point to note is that the values of seem to be slightly 
larger than would be expected from the AFM imaging, 
which is likely connected with the fact that the AFM tip 
has a finite size and cannot go into grooves]. 

We thus conclude that the description of flow near 
rough surfaces has to be corrected, but this is not the 
relaxation of no-slip boundary conditions, but the cor- 
rection for separation with its shift to larger distances 
within the asperities size. It has to be stressed that sim- 
ilar ideas were justified theoretically at macroscale for 
a shear flow along periodically corrugated wall [13 and 
recently for a far-field motion of a sphere towards such 
a wall [2^. They are also consistent with molecular- 
dynamic simulations on simple model systems Here 
for the first time we provided accurate experimental data 
showing that the concept can be applied for a randomly 
rough wall and at the nanoscale, down to a contact. Note 
that Eq.[7|has not received so far any theoretical justifica- 
tion for short, i.e. of the order of or smaller than the size 
of asperities, separations. This is probably because this 
situation escapes from the framework of the lubrication 
approximation since two length scales of the problem be- 
come comparable. The applicability of Eq.[7|at the short 
distances probably reflects the fact that the height of the 
roughness elements is much smaller than the sphere ra- 
dius, so that even when separation is getting small, the 
plane "sees" many roughness elements at the same time 
and fluctuations are averaged out. This hypothesis re- 
mains the subject for further theoretical work. 

This conclusion remains valid for all driving speeds, 
and results are summarized in Fig. One can see that 
rs remains roughly the same at all speeds, although it 
shows some weak tendency to increase with the rate of 
approach. However, this increase is within the error of 
experimental data, so it remains an open question for 
future research. The same remarks are true for bg, which 
is included for completeness since slip-model is applicable 
at large separations. 

Note that our results clarify the reason for appar- 
ently contradictory reports suggesting that roughness in- 
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FIG. 5: A hydrodynamic height of asperities hs (open cir- 
cles) and a hydrodynamic shift Ts (filled squares) obtained at 
different speed of approach. 



crease 0| and decrease ^3 the drag force. In our opin- 
ion these different conclusions only reflect that in these 
papers the wall location was defined differently. It fol- 
lows from our results that if equivalent smooth surface is 
defined as coinciding with the location of valleys of as- 
perities ^hich corresponds to experimental techniques 
used in ,14]), the measured force is larger than expected 
for an equivalent surface (cf. dotted curve in Fig. ^ . If, 
however, it is defined on the peaks of asperities (experi- 
mental techniques used in [l7|), then the measured force 
is smaller (cf. dashed curve in Fig. Clearly, both 
results [ij, [13 physically mean that roughness increase 
the dissipation in the system, and that an equivalent sur- 
face is located at the intermediate position. Note, that 
in {12]. which reports the increase in force with rough- 
ness, the equivalent surface is also defined on the peaks 
of asperities. We do not have any explanation of this 
result. 

Our data and conclusions do not confirm results ob- 
tained with similar, but "smooth" surfaces, where shear- 
rate dependent slippage was detected |2j]. In our 
opinion, the reduction in hydrodynamic resistance force 
might indicate that their gold-coated thiolated sphere 
was in reality rougher than expected. Reasons for 
a shear-dependency could probably be connected with 
some errors in experimental determination of dh/dt since 
the piezotranslator used in and later in 0|, is not 
suitable for highly dynamic force measurements due to 
its non-linearity. Another reason for a significant rate- 
dependence is the use of binary mixtures (U]. Clearly, 
both confinement and shear might lead to their stratifi- 
cation and a formation of a thin low viscosity lubricant 
layer |inj |. It is well known and has been proven exper- 
imentally that such a layer leads to rate-dependent 
phenomena (normally expected only at a very large shear 
rate [2^) even at low speed. We suspect that this effect 
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of stratification of binary mixture, enhanced by rough- 
ness, is responsible for a very large reduction in force 
observed in Finally, we would like to stress that 

since the force balance represents a differential equation, 
even small could implicate the results by decreasing 
dh/dt, and therefore Fh- The similar remark concern Fg, 
which cannot be excluded from analysis. Both Fg and Fc 
are not present in the force balance specified in 0, |^ , 
which might cause further inaccuracy in results. 

In summary, by performing high-speed drainage exper- 
iment with nanorough hydrophilic surfaces we demon- 
strate that their interaction is similar to that between 
equivalent smooth surfaces located at some position be- 
tween peaks and valleys of asperities. Thus, our results 
are in favor of no-slip assumptions for a hydrophilic sur- 
face valid down to a contact. 
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